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Using (infinite) density matrix renormalization group techniques, ground state properties of anti- 
ferromagnetic S=l Heisenberg spin chains with exchange and single-site anisotropies in an external 
field are studied. The phase diagram is known to display a plenitude of interesting phases. We 
elucidate quantum phase transitions between the supersolid and spin-liquid as well as the spin- 
liquid and the ferromagnetic phases. Analyzing spin correlation functions in the spin-liquid phase, 
commensurate and (two distinct) incommensurate regions are identified. 

PACS numbers; 75.10.Jm, 75.40.Mg, 75.40.Cx 



In recent years, ground state properties of the antifer- 
romagnetic Heisenberg spin~l chain with single-site and 
uniaxial exchange anisotropies in an external field have 
been investigated rather extensivelj^"— . The model is de- 
scribed by the Hamiltonian 
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where i denotes the lattice sites, A the exchange, and 
D the single-site anisotropy. The external field B acts 
along the ^-direction. 

The magnet displays various intriguing phases at 
zero temperature (and, thence, corresponding quantum 
phase transitions), the antiferromagnetic (AF), ferro- 
magnetic (F), half-magnetization plateau (HMP), spin- 
liquid (SL), supersolid (SS), and 'large-D' phases. Some 
of these phases, the AF, F, SL, and SS phases, show up in 
the corresponding classical Heisenberg model'^'^'S, while 
the HMP and large-D phases reflect the discretization of 
the spin orientations in the quantum case. The theoret- 
ical efforts have been motivated and inspired, partly, by 
related experiments^'*. 

Perhaps most attention, in the context of this, Eq. (1), 
and similar—"— models, has been paid to the supersolid 
phasei^, being the analog of the 'mixed' or 'biconical'— 
phase in the classical limi t^^'^* . Note that a mapping 
from quantum lattice gases to magnetic systems, explain- 
ing the term 'supersolid' for magnets, has been intro- 
duced some decades agoi2. Typically, quantum fluctua- 
tions tend to reduce substantially the range of stability 
of the supersolid phase^!^, as compared to the classical 
variant. 

The spin correlations in the supersolid phase of the 
anisotropic Heisenberg spin chain, (1), have been argued 
to behave like in a Luttinger liquid, with algebraic spatial 
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FIG. 1. Staggered magnetization nist vs. external field B/J 
at D/ J = A/2 — 2.5 near the supersolid to spin-liquid tran- 
sition. For comparison, a power-law fit to the IDMRG data 
with the critical exponent /? = 1/8 is shown, see text. 



decayS,. Magnetization proflles have revealed the close 
analogy of the supersolid to the corresponding classical 
biconical structures^i^. The critical exponent of the spin 
stiffness, describing the transition to the bordering AF 
and HMP phases, has been found^ii^ to be 1/2. 

In this contribution, we shall consider interesting as- 
pects of the model (1) which have not been studied in de- 
tail so far. We shall deal with the transition between the 
supersolid and the spin-liquid phases as well as with the 
SL-F quantum phase transition. In addition, spin corre- 
lations in the SL phase will be analyzed, to clarify, espe- 
cially, previous suggestions on distinguishing commensu- 
rate and incommensurate regions in that phasei*^"^. 

In the present study, mainly infinite density matrix 
renormalization group (iDMRG) techniques?^"?? have 
been used, with systematic enlargening on the number 
of matrices in the matrix product states. In general, the 
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FIG. 2. Longitudinal spin correlation function [Fz — rri^], 
vs separation distance r, at A = 5, with (a) D/J — 2.5, 
m = 23/40 (circles), (IC)i type, (b) D/J = -1.5, m = 1/5 
(squares), (IC)i type, and (c) D/J = 1, m = 4/10 (dia- 
monds), C type. Interpolating lines are guides to the eye. 



chosen size of the matrices is an important parameter de- 
termining the rehability of the calculation. Here, the di- 
mension M of the largest matrices ranges, typically, from 
50 to 500. The truncation error varies in between 10~® 
and 10~^°. In a few cases, results are compared to ones 
we obtain from DMRG calculations for finite chains, with 
open boundary conditions, of length L, with L < 128. 

Following previous analysea^"— of the model (1), we fo- 
cus on two cases: at fixed ratio between the two types 
of anisotropics, D/J = A/2, and at given, quite large 
exchange anisotropy, A = 5, with varying single-site 
anisotropy D/J. (Parts of) the corresponding ground 
state phase diagrams have been obtained before^"— , us- 
ing DMRG and quantum Monte Carlo techniques. In 
both cases, the phase diagrams include the AF, F, HMP, 
SL, and SS phases. The supersolid phase results from 
competing uniaxial, along the z-axis, exchange, A > 1, 
and planar single-site, D > 0, anisotropics. 

The supersolid phase may be bordered by massive, AF 
or HMP, or by critical, SL, phases^Ti^. The transitions 
to the massive phases have been investigated in detail 
before^.!^. Here we discuss the transition from the su- 
persolid to the spin-liquid phase, the SS-SL transition. 
As illustrated in Fig. 1 for A = 2D /J = 5, the tran- 
sition seems to belong to the two-dimensional classical 
Ising universality class. The critical exponent /3, describ- 
ing the vanishing of the staggered magnetization mst at 
the transition is, indeed, consistent with the famous On- 
sager value /3 = 1/8. Actually, we did a x^-fit of our 
iDMRG data in the range 7.217 < B/J < 7.262 to the 
form rUst = a{Bc/J - B/Jf/^. We obtain B^ « 7.261, 
as shown in Fig. 1. Note the (small) deviation extremely 
close to Be, which one may attribute to discretization 
error in determining the magnetic field or to insufficient 
size of the matrix dimension M, M < 500, in this regime 
near the transition. 

In addition, we also identify an Ising-like sector near 



the SS-SL transition, with exponentially decaying lon- 
gitudinal spin correlations up to rather large distances, 
as will be discussed below. Note that our suggestion 
on the universality class of the SS-SL transition is in 
line with a corresponding finding on a related two- 
dimensional quantum anisotropic Heisenberg antifcrro- 
magnet, where the supersolid to spin-liquid transition 
has been concluded to be in the universality class of the 
three-dimensional Ising case^i. It is worth mentioning 
that both suggestions, for quantum magnets in dimen- 
sions d= 1 and 2 at zero temperature, agree with the 
well-known dimensional argument^ relating critical ex- 
ponents in d-dimensional quantum systems to those in 
corresponding (d + l)-dimensional classical systems. In- 
deed, here the classical transitions between the biconical 
and spin-flop phases are of Ising type^ii^. 
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FIG. 3. Solid boundary lines separate C and IC regions in 
the {D/J,m) plane at A = 5. The dashed line sketches, as 
a guide to the eye, the border of the forbidden region (FR), 
see Ref. 1. The dashed-dotted line in the (IC)i phase di- 
vides ferroquadrupolar, at high m, from spin-density-wave 
ordering. 

Before turning to the discussion of the SL-F quantum 
phase transition, let us first consider characteristics of 
the SL phase. As has been noted before^'^'^, the Hamil- 
tonian (1) may describe both commensurate (C) and in- 
commensurate (IC) spin-liquid structures, as has been 
inferred from the behavior of energy gaps^ and magneti- 
zation profiles^i^. 

In this study, we shall present direct evidence for both 
types of structures in the spin-liquid phase by analyzing, 
especially, longitudinal Tz{r) = (S'fS'f^.^) spin correla- 
tion functions. Asymptotically, r — > oo, F^ acquires 
the value m^, where m is the total magnetization per 
site. For sufficiently large distances, r, the dominant 
algebraically decaying term of the correlations may be 
expected^* to be of commensurate (C) form, cx 1/r^, 
or of incommensurate (IC) form, oc (1/r'') cos(gr), with 
T] < 2, as usual for Luttinger liquids. Such a behavior is 
confirmed by our iDMRG calculations. In the IC case, 
we find the wavenumber q to be related to the total mag- 
netization per site, m, in two distinct ways: We obtain 
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either qi = 7r(l - to) ((IC)i), or ga = 27r(l - m) ((10)2), 
setting the lattice constant equal to one. 

Examples of longitudinal correlation functions of type 
C, (IC)i, and (IC)2 are depicted in Fig. 2, for selected 
values of D/J and to, fixing the exchange anisotropy, 
A = 5. Note that in the example for the C case, 
shows roughly an exponential decay with oscillations at 
small distances r, approaching the monotonic algebraic 
decay, oc only at larger separations. 

Varying systematically the single-site anisotropy D 
and the magnetization to, at A = 5 (compare to Ref. 
1), one may then identify three different regions, C, 
(IC)i, and (10)2, in the {D/J,m) plane, as shown in 
Fig. 3. At sufficiently small negative values of D, one 
observes the (IC)i region. This region may be subdi- 
vided into two parts: At larger magnetizations, we ob- 
serve ferroquadrupolar orderingSS^-!^^ ^ where the algebraic 
decay of the four-point transverse correlation function 
{{StyiS~+r)'^) is slower than that of F^(r), due to a 
smaller exponent rj. At lower magnetizations, one en- 
counters a spin-density-wave ordering, with the longitu- 
dinal spin correlations being dominant. The (IC)2 region 
occurs at larger values of D/J. In between the two IC 
regions, the commensurate region intervenes. There, the 
exponent rj characterizing the algebraic decay of with 
spatially modulated behavior, is larger than 2. Asymp- 
totically, for large distances r, the dominant algebraic 
term is then proportional to decaying monotoni- 

cally. Indeed, the changeover between the C and IC re- 
gions may be conveniently monitored by determining the 
exponent rj from fits of the iDMRG data for the longi- 
tudinal spin correlations^^. A few examples, at several 
fixed values of D/J and changing to, are displayed in Fig. 
4, compare to Fig. 3. 

The 'forbidden region' (FR) in the {D/J,m) plane, 
which has been sketched in Fig. 3, gives rise to first 
order transitions. 

Our calculations on the spin correlations in the 
(Z?/J, to) plane confirm and refine substantially previ- 
ous findings^, where the (IC)2 region seems to have been 
overlooked. Similarly, there has been no mentioning of 
the two distinct parts of the (IC)i phase and of the super- 
solid phase which shows up at fairly low magnetizations 
and sufficiently large values oi D/J, as we have discussed 
before^. We omitted the supersolid phase in Fig. 3, for 
reasons of simplicity. 

Close to the supersolid to spin-liquid transition, 
the longitudinal spin correlations are governed, up to 
quite large distances, by an exponential decay with tt- 
oscillations, signalling, presumably, the above mentioned 
Ising sector. Further details are presented elsewherSiSZ . 

We now turn to the discusssion of the SL-F transition. 
The phase transition may be characterized by the behav- 
ior of the total magnetization per site to, as illustrated 
in Fig. 5 for Z? = 0, with the SL phase being of (IC)2 
type. Obviously, the (i-)DMRG data may well be fitted 
to the form 
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FIG. 4. Exponent i]^ for the oscillatory part of the alge- 
braically decaying longitudinal spin correlations vs. mag- 
netization m at A = 5 and various values of D/J. 
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FIG. 5. Magnetization m versus field B close to the SL-F 
transition at A = 5 and D = 0. Data from iDMRG and 
finite-size DMRG, for L— 128 sites, calculations are shown, 
together with a square-root power law fit (solid line), see (2). 



l-moi{Bc/J - B/jy/^ (2) 

where Be is the critical field of the SL-F transition. 
Indeed, we did a x^-fit of the (i)DMRG data in the 
range between 11.7 < B/J < 12.0 to Eq. (2), deter- 
mining the proportionality factor and the critical field 
Be ~ 12.0. Further away from the transition, devia- 
tions from the simple power law may be observed, see 
Fig. 5. Note that in the (IC)2 region of the SL phase, 
one has 1 — to oc g. Thence, Eq. (2) corresponds to 
the well known Pokrovsky-Talapov^* square-root power 
law for the wavenumber q, describing the C-IC tran- 
sition in two-dimensional classical systems with uniax- 
ial spatial anisotropy. Indeed, it seems tempting and 
reasonableSSiSS to associate the SL-F transition for the 
quantum spin chain with that universality class. 

In summary, studying ground state properties of the 
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S=l anisotropic Heisenberg antiferromagnetic chain, us- 
ing mainly the iDMRG approach, we present evidence 
for the spin-hquid to supersohd transition being in the 
two-dimensional Ising and for the SL-F transition being 
in the Pokrovsky-Talapov universality classes. The spin- 
liquid phase is found to consist of commensurate and two 



distinct incommensurate regions. 
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